Fusion rules and vortices in p x + ip y superconductors 
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Abstract 

The "half-quantum" vortices (er) and quasiparticles (tp) in a two-dimensional p x + ip y supercon- 
ductor obey the Ising-like fusion rules i/jxip = l,axip = a, and a x a = 1 + ip. We explain how the 
physical fusion of vortex-antivortex pairs allows us to use these rules to read out the information 
encoded in the topologically protected space of degenerate ground states. We comment on the 
potential applicability of this fact to quantum computation. 



PACS numbers: 71.10.Pm, 73.43.-f, 74.90.+n 



I. INTRODUCTION 

The magic of a quantum computer jjj, |2j, |3|, |4| is that it makes a virtue out of a vice. The 
numerical simulation of even simple quantum systems requires exponentially large storage 
space and exponentially long computation times. The quantum state of a collection of iV 
spin-i particles, or qubits, requires 2 N classical variables for its description, and solving 
the Schrodinger equation to follow its time evolution requires a conventional computer to 
perform up to 2 N x 2 N operations per time step. Conversely, a computer built out of 
quantum components would be able to exploit the massive parallelism inherent in quantum 
time evolution to provide fast solutions for problems that would require exponential time 
on conventional machines. 

The difficulties to be overcome in building a quantum computer are many. The quantum 
system at its core must be strongly isolated from the environment so as to avoid decoherence; 
the unitary transformations that constitute the elementary computational steps must be 
performed with sufficient precision that error-correcting codes p| remain effective; and, 
after all the quantum computation has been performed, the system must be capable of 
being reconnected to the outside world in such a way that the output can be read off via a 
measurement process that does not disturb the result. 

A particularly appealing scheme for quantum computation exploits topologically pro- 
tected macroscopic quantum states in many-body systems 0, Q, 0, S| . The first examples of 
such topologically ordered and protected many-body quantum states were found by Xiao- 
Gang Wen and Qian Niu n| in the context of the fractional quantum Hall effect. Wen and 
Niu observed that if we place a filling-fraction u — 1/n quantum Hall state on a genus g 
Riemann surface (a thought experiment only!) then there are n 9 degenerate, and essentially 
indistinguishable, ground states. They also observed that by creating a vortex antivortex 
pair in the Hall fluid and then moving the vortex around one of the homology generators of 
the surface they could cause the system to roll over from one degenerate state to the next. 
Since the degeneracy is lifted only by the exponentially suppressed tunneling of vortices 
around the generators, the coherence of a linear superposition of such states is topologi- 
cally protected. Furthermore, as described above, such superpositions can in principle be 
manipulated by controlling the motion of their vortex defects, and so inducing on them 
representations of the braid group of vortex world-lines. 
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For simple systems, such as the v = 1/n Hall state on a simply-connected surface, these 
braid-group representations are abelian, interchange and braiding giving rise only to a phase 
factor. This is not sufficient for a universal quantum computer ]8Llll|. We need a non- 
abelian representation — i.e. particles with non-abelian statistics [l2j. There are a various 
two-dimensional physical systems that should, in theory, have excitations with non-abelian 
statistics. Those closest to being realized in practice are all associated with many-body 
wavefunctions that contain a Pfaffian factor. These candidates are: i) a quantum Hall effect 
state seen at filling-fraction v = 5/2 ii) an anticipated, but as yet unobserved, phase of 
a rapidly rotating atomic Bose gas, iii) a two-dimensional p x + ip y superconducting state— 
an example of which has possibly been observed in the layered Sr 2 Ru04 superconductor 
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151 ] - While the non-abelian statistics of these Pfaffian states is still not quite adequate 



for the construction of a universal computer, it is still worthy of study. 

The mathematical origin of the non-abelian braiding in the first two candidate systems is 
quite technical, involving complicated linear dependences among various Pfaffian wavefunc- 
tions ^| . The non-abelian braid statistics of the p x + ip y superconductor is much easier to 
understand physically. The mechanism has been beautifully explained by Ivanov and 
further clarified by Stern et al. Q|. 

The essential ingredient is that the core of an Abrikosov vortex in any superconductor 
contains localized low-energy bound states whose existence is guaranteed by an index theo- 
rem that relates the phase winding number n of the vortex to the number of branches of low 
energy excitation s |19^ . The p x + ip y superconductor is unusual in that it can host a "half- 
quantum" vortex |2Cj, where the phase winding is seen by only one of the two components of 
the electron spin. The quasiparticles bound in the core of a two-dimensional version of such 
a vortex are Majorana, i.e. they are identical to their antiparticles. Further, for odd phase 
winding numbers, each vortex binds a zero energy mode whose field-expansion coefficient 
is hermitian. In the presence of 2N such vortices, the many-body Hilbert space of zero- 
energy states is 2 7V_1 dimensional. This number is in contrast with the naive count of 2 2N 
dimensions that would come from (wrongly) supposing that each single-particle zero-energy 
mode must be either occupied or unoccupied. Because these 2 N ~ 1 independent states can- 
not be associated with any individual vortex, the many-body state encodes highly non-local 



information. Ivanov 17( demonstrated that the operation of braiding or interchanging the 
vortices acts on this zero-energy space in the same way that braiding acts on the 2 Ar_1 di- 
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mensional space of iV-point conformal blocks in the level k = 2, SU(2) Wess-Zumino-Witten 
(WZW) model. 

It is not surprising that this particular WZW model plays some role. It is known that it 
describes the low energy physics of Pfaffian quantum Hall state |2l|. It also describes the 

23| . and is a component of the low 



low energy physics of the Bose gas pfaffian state 

energy effective action of the two-dimensional p x + ip y superconductor 24] . 

Because the zero-energy states consist of an equal superposition of particle and hole, 
they are electrically neutral. They are also nonlocal. As a consequence they should be 
little affected by local impurities and external fields, and so be topologically protected. 
The vortices binding them, however, carry magnetic flux and so might manipulable by 
external probes such as STM tips. These vortices are thus potential candidates for the 
basic building blocks of a quantum computer. The problem is that the internal state in the 
2 iV ~ 1 -dimensional protected zero-energy space is so decoupled from the outside world that 
it is hard to see how it can be measured. 

This read-out problem is not unique to non-abelian braid statistics. There is as yet no 
completely convincing experimental demonstration of even the abelian "anyonic" statistics 
in the v — 1/n Hall state — although no theorist doubts that it is present. An experiment 
that can detect non-abelian statistics is even harder to perform. We need some method of 
accessing the non-abelian state by reading out the information in the protected space. 

One approach, at least in principle, is to use the "fusion rules" of the vortices to re- 
connect the protected states to something measurable. The vortex fusion algebra contains, 
in addition to the obvious rule that phase- winding number adds, a Z2 factor that is essentially 
that of the Ising model primary-field operators I, a, and ip [25], or of the analogous affine 
primary fields (of spin zero, one-half, and one respectively) in the SU(2)2 WZW model: 

■0 X if) = I, a x a = I + ip, fxff = a. (1) 

In the superconductor it is natural to identify I with the ground state and ip with the 
Bogoliubov quasiparticle. The latter is conserved only mod 2 because two quasiparticles 
can combine to form a Cooper pair and vanish into the condensate. Hence ip X ip — > I. 
Slightly less obvious is the identification of o with the odd-winding-number vortex. This 
comes about because a pair of vortices is associated with a single zero-energy state that can 
be occupied (?/>) or unoccupied (I). Hence a x a = I + ip. The third product, a x ip = a, 



is then determined by the associativity of the fusion-product algebra. The principal claim 
of this paper is that it is possible to use this algebra to determine some of the information 
in the topologically protected space by physically fusing a vortex with an antivortex and 
seeing if they leave behind a real quasiparticle ip. 

In subsequent sections, we will fill in the details of this process. In section two we 
discuss the form of the zero energy solutions to the single-particle Bogoliubov-de-Gennes 
equation, focusing on the sign ambiguities which are the ultimate source of the particle- 
creation legerdemain. In section three we show how tunneling between nearby vortices can 
be used to resolve the sign ambiguity, and how braiding followed by fusing reveals details 
of the protected state. Section four is a conclusion and discussion. An appendix describes 
some necessary, but technical, aspects of Berry transport in a superconducting state. 

II. VORTEX CORE STATES 

The topological objects that have the properties we seek to exploit are the "half-quantum" 
vortices of a thin film of p x + ip y spin-triplet superfluid, such as 3 He-A, or possibly the 
superconductor S^RuO^ The order parameter in such a superfluid contains an angular- 
moment vector 1, a vector d characterizing the Cooper-pair spin state, and an overall phase- 
factor exp(zx). In a half-quantum vortex the vector 1 lies along the vortex axis and the 
vector d lies in plane which, for convenience, we imagine to be perpendicular to 1. As we 
encircle the vortex the phase x increases by 7r while the d vector rotates through 180°. 
Although exp(zx) and d separately change sign, it is their product that appears in the order 
parameter, and this is single-valued. The effect of the combined rotation is that the spin- 
up fermions see the order parameter phase wind through 2ir, while the spin-down fermions 
see no phase change at all. Andreev reflection off the winding phase will bind low energy 
spin-up quasiparticles modes in the vortex core, but there will be no spin-down quasiparticle 
states with energy significantly less that the gap |A|. Because only one component of spin 
is involved in the low energy physics, we will from now on regard the fermions as being 
spinless. It is this spinlessness that allows the p x + ip y superconductor to have topological 
properties analogous to those of the spin polarized Pfaffian quantum Hall state . 

Consider a two-dimensional film of superfluid with its 1 vector perpendicular to the film, 
and with a winding-number n vortex at the origin. We are principally interested in the 
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case n = ±1, but we will keep n general for the moment, so as to bring out a distinction 
between vortices with n even and n odd. The Bogoliubov-de-Gennes (BdG) equation 
has bound-state solutions of the form 




e™[ Ul ) (2) 



where 




(3) 



e ie(n+1 V 2 [a(r)H l+1/2 (k f r) + c.c] 
e- i6( - n+1 ^ 2 [b{r)H^ 1/2 {k f r) + c.c] 
Here Hi(kfr) is a Hankel function of the radial co-ordinate r, the angle 9 is the polar co- 
ordinate, and kf is the Fermi momentum. The coefficient functions a(r), b(r) are slowly 
varying on the length scale of kj 1 and are most conveniently written in terms of a pair of 
auxiliary variables x and 8(x). These are defined in terms of the impact parameter r = l/kf 
by 



x = \ r 2 — Tq, (4) 



6{x) = tan" 1 ( y ) . (5) 



(6) 



If we write 

a(r) \ ( e- ind ^/ 2 a(x) 
b(r) ) \ e ine{x)/2 b{x 
then, in the Andreev approximation, the coefficients a(x) and b(x) obey 

-iv f d x A(r)e in9 ^\ f a(x)\ ( a(x)\ 

A(r)e- in6 W +iv f d x ) \ b(x) / ~ \ b(x) ) ' 

The one-dimensional eigenvalue equation ((7j) has a physically appealing interpretation as 
describing the propagation of a quasiparticle along a rectilinear trajectory with ro its distance 
of closest approach to the vortex centre. The variable x is the distance along the trajectory 
with x = being the point of closest approach. As the quasiparticle moves it sees the 
local value of the order parameter A(r(x)) exp[in6(x)]. The Hankel function is evanescent 
for r < r , and a and b can be taken as being constant in this region. In order for the 
Hankel functions to combine to give the Bessel functions Ji±i/2{kr) that remain finite at 
r = 0, we need to impose the condition that a and b (and hence a and b) be real at x = 0. 
This condition allows us to extend the definition of a(x) and b(x) to continuous functions 
on the entire real line by setting a(— x) = (a(x))* and b(— x) = (b(xj)*. It may be verified 



that (a(x),b(x)) T continues to satisfy (JZJ) in the extended domain. The one-dimensional 
equation is able to capture the physics of the bound states because the Andreev scattering 
that confines the particle/hole in the vortex core is almost exactly retro-reflective J^. 

Solving the one- dimensional problem shows that, for small /, the energy eigenvalue is 
given by E(l) = —Iujq where ujq is a frequency determined by the order parameter profile 
A(r). This frequency is positive if n > and negative if n < 0. A glance at (J2J) and (JHJ) 
shows that single- valuedness of the BdG wavefunction requires / to take integer values when 
n is odd, and half-integer values when n is even. The spectrum may therefore be labeled by 
an integer m and is given by 

{—u m, n odd 

( 8 ) 
— cu (m + 1/2), n even. 

Both cases are consistent with the E «-> — E symmetry of the BdG eigenvalues. For odd 
winding numbers the spectrum contains an E = mode. If we create N vortices there will 
be one of these zero modes localized in each vortex core. There will, however, always be an 
even number of zero modes. For an odd number of vortices, an additional zero-mode will 
be found on the boundary of the superfluid 

Although the zero modes for the n = +1 and n = — 1 vortices both have / = 0, their 
actual angular dependence differs. For the n — +1 vortex, whose circulation is in the same 
sense as the p x + ip y Cooper-pair angular momentum, we have 

u (r,9) 



vo(r,6) 




(9) 



For the n = — 1 vortex we have, 





(10) 

In both cases the radial function /(r) is real and decays as exp(— Ar/vf) away from the 
vortex core. 

We have been tacitly assuming that A is real. A global phase rotation A — > e iX A alters 
each of the eigenmodes and we can choose an overall phase for the modes so that 

u l {r,6)\ ( e ix / 2 ui(r,9) 



vi{r,6) J \e- ix /\(r,6) 



(11) 
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An inspection of the detailed form of the solutions shows that, with this phase choice, we 
retain the property that the upper component of eigenmode I is the complex conjugate of 
the lower component for eigenmode — /. This property does not uniquely specify the {u,v) T 
vectors, however. Even after normalization, so that J (\ui\ 2 + \vi\ 2 ) cPx = 1, a choice of overall 
sign remains to be made. 

In the vicinity of the vortex core we can use the bound states to make a mode expansion 
of the low energy part of the fermion field 

$(r,0)\ ^ fu t \ 

J — y b\ J e l + higher-energy modes. (12) 
y{r,B)J , \vij 

The fact that the lower field component is the hermitian conjugate of the upper compo- 
nent ^ coupled with the phase choice made in (jllj) . enforces the hermiticity condition bi = b_ t 
on the annihilation and creation operator coefficients. This property is characteristic of a 
Majorana fermion. The bound-state quasiparticles are therefore their own antiparticles. For 
n odd, the zero-energy, mode coefficient obeys 60 — ^0 an d cannot be thought of as being 
either a creation or annihilation operator. Instead, the conditions {^(x), ip(y)} = and 
{$(x), ip(y)} = 5 2 (x — y) together with the completeness of the eigenmodes tell us that 

6g = l/2, {6o,6,} = {6o,6j} = 0, 1^0. (13) 

We have taken the trouble in this section to display the bound-state wavefunctions in 
some detail. We did this to make clear the origin of the b\ = b-i Majorana condition, 
to make explicit our phase choices, and to point out the remaining sign ambiguity. This 
last point is important. For conventional fermion annihilation and creation operators, the 
relative sign of a« and a\ is fixed by the condition {a,a^} = 1, but the relative sign of the 
various eij with respect to one another is unimportant as it can always be absorbed into the 
definition of the free-field modes whose coefficients they are. The relative sign of the b 's 
associated with different vortices has a physical significance that will be appreciated after 
we consider tunneling between nearby vortex pairs. 



III. BRAIDING AND FUSION 



Suppose we start from a homogeneous (A = const.) state and adiabatically create iV 
vortex-antivortex pairs. For ease of discussion imagine the the pairs arranged so that each 
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antivortex lies on the x axis, and the corresponding vortex lies vertically above. We label 
the pairs from left to right by the index i — 1, . . . , N. In the core of each vortex and 
antivortex there will be a zero-mode with its corresponding b mode operator. Up to a 
factor, these hermitian operators obey gamma-matrix anti-commutation relations, so it is 
natural to rename these 6 ' s as li/V^ (vortex) and and ji+iy/V^ (antivortex). These 2N 
operators can be assembled to make iV each of annihilation and creation operators 

h = ^(7. + ili+N), 

A = \{li-ili+N)- (14) 

These fej and fej's act on a 2 Ar -dimensional space of degenerate ground states. This space 
is split into two physically equivalent 2 Af ~ 1 -dimensional spaces by a superselection rule: 
although fermion number conservation is broken by the presence of the condensate, no 
physical process internal to the system can change an odd total fermion number into an 
even fermion number. Which of the two spaces we are in depends on whether the total 
number of particles in the system is odd or even. 

The odd-even decomposition is reflected in the fusion algebra. By using the rules (JTJ) to 
repeatedly fuse a o with a x a = I + ip, we find 

a x a x a = a + a, 
a x a x a x a = I + ip + 1 + ip, 

etc. In general, for 2N vortices, 

a x ■ ■ ■ x a = ( I + ■ ■ ■ + I ) + ( i/j + ■ ■ ■ + t/j ) . (15) 

2 N factors 2 N ~ 1 terms 2 N ~ 1 terms 

The 2 7V ~ 1 -dimensional multiplicity spaces of the I and the ip are separately invariant under 
the action of the braid group. 

When each vortex "i" is brought close to its sibling antivortex U N + the 2 Ar -fold 
ground state degeneracy is lifted. Because two nearby odd-winding number vortices are 
effectively a vortex with an even winding number, tunneling between their two zero-energy 
BdG modes will cause them to realign and become two eigenmodes with some small non- 
zero energy ±£?j. In the second quantized Bogoliubov Hamiltonian these two single-particle 
modes make one state that can either be occupied and contribute to the total energy 
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or be unoccupied and contribute energy —Ei. If we were to slowly merge each vortex with 
its sibling antivortex, the ±Ei eigenmodes will split further and eventually merge with the 
continuum of unbound states. If the Ei state is empty, nothing remains after the vortices 
annihilate. If it is occupied, an unbound quasiparticle with energy of |A| is left behind. 

It is not easy to compute the matrix elements for the tunnel splitting, but we do not 
need to know them in detail. When the splitting is small we can restrict ourselves to the 2 N 
dimensional space of nearly degenerate states, and in this space the many-body Hamiltonian 
must be of the form 

N 1 1 N 

^ = Z^4-[7JV+i,7i] + Ei 
i=i i=i 

N N 

i=l i=l 

N 

= (16) 

i=l 

This is because -^[jN+i, 1i] is the only hermitian operator that can be made solely out of 7, 
and 7i+iv, and so the only possible operator that tunneling between vortex "z" and U N + z" 
can contribute to the Hamiltonian. Hermiticity demands that that the Ei be real. We 
can also assume that the Ei are positive by using this requirement to fix the relative sign 
ambiguity between •ji+N and 7$. 

The lowest energy state |0) for an even number of particles is that which is annihilated by 
all the hi. In this ground state all fermions are paired, and all of the not-quite- zero-energy 
states are empty. If the system contains an odd number of particles, however, one will be 
left unpaired, and will occupy the lowest of the not-quite-zero-energy states. 

Now we consider how we can manipulate the occupation numbers of the not-quite- 
zero modes. Following Ivanov we adiabatically transport vortex % around vortex 
j and bring it back to its original position. In this process the local phase \ seen by 
each vortex will increment by ±27r, and so cause the phase factors e* x ^ 2 (see (llljl ) in 
the zero modes of vortex i and j to change sign. The field operators ip(x), tp^(x) are 
indifferent to the choice of modes in which we expand them, and must be unchanged 
by the braiding process. The sign change of the mode vector (v,q, Vq) must therefore be 
compensated by a change in the sign of the mode coefficients 7* and 7,-. Consequently 
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FIG. 1: Phase changes due to the braid operation T±: Before the braiding the phase X2(l) seen 
by vortex 1 due to vortex 2 is just more than 180°, and the phase Xi(2) seen by vortex 2 due to 
vortex 1 is just more than 0°. After the braiding the phase Xi(2) at vortex 2 due to vortex 1 is 
just more than 180°, but because vortex 1 remains below below the displaced vortex-2 branch cut, 
the phase X2(l) is now just more than 360°. Consequently 71 is replaced by 72, but 72 is replaced 
by -71 • 



Ho = --- + E~[ lN+i , 7i ] + • • • + Ej-[y N+j , 7i ] + • • • + + • • • + -Ej + ■ 

is changed to 

= . . . _ E~[ lN+i , 7i ] + ■ ■ • - E~[ lN+j , 7i ] + • • • + lEk + • • • + \Ej + 



it acquire a non-abelian, 
~l 

18| , that, for the choice 



In addition to this explicit monodromy in Hq, the state |0) mig 
holonomy from Berry transport. It was argued by Stern et al 
of phases in (u ,v ), the state |0) is at most multiplied by an overall Berry phase, and so 
essentially returns to itself after the braiding process. We agree with this conclusion, al- 
though we find the discussion in unnecessarily involved. We therefore provide our own 
derivation of this key fact in the appendix. 

After the braiding the state |0) remains an eigenvector of all the j^lN+kilk] with eigen- 
value — 1/2. Because of the sign changes in H Q , however, it is no longer the lowest energy 
state of Hq™. It is instead an excited state with energy E t + Ej, corresponding to the 
two not-quite-zero modes of the i and j-th pairs being occupied. If we now adiabatically 
fuse the i-th vortex with its sibling antivortex we recover the uniform state together with a 
quasiparticle ip. The same is true for for the j-th pair. 

We now suppose we make an anti-clockwise interchange Tj of the i and % + 1-th vortices. 
In order to follow what happens, we must first fix the sign ambiguity between the 7i . We 
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therefore draw "branch cuts" to the right of each vortex and parallel to the x-axis and set 
the Xi produced by the vortex equal to zero immediately above the cut. The local phase 
x{i) seen by vortex % is then the sum of the Ylj^i Xj' s °f the other vortices and anti-vortices. 
Under the interchange, and keeping track of the explicit monodromy of the local x«' s ) we 
find that the zero mode of the % + 1-th vortex replaces that of the i-th while the z'-th zero 
mode replaces minus that of the z + 1-th. Thus (see fig ^ and 3|) 

r ii -> n+i , x 

Ti : (17) 
I li+i -> ~li- 

Since the geometric arrangement of the vortices is unchanged, the Hamiltonian becomes 

F new = ■ • • + E~[ lN+i , 7i+1 ] + E i+1 ^[j N+i+1 , 7i ] + ■ ■ ■ + X -E { + • • ■ + ^E i+1 + ■■■. 

The crucial effect is not so much the sign change but that the Berry transport of |0) preserves 
its property that it is the state killed by the bk of the original r=[lN+k,lk]- Consequently 
|0) is no longer an eigenstate of the new Hamiltonian, but is instead a linear superposition 
of eigenstates. The outcome of fusing vortex z with its antivortex is no longer certain. Same 
is true for the outcome of fusing vortex z + 1 with its antivortex. In fact, it will be shown 
later that we have constructed an entangled state. 

The sign change in 7^ —> — jj does have significance, however, as it ensures that the result 
of braiding is to take 7, —> 7^ with 

li = iPji = UnU-\ 

where O^- is an element of SO(2iV), as opposed to 0(2N). Here U is a spin-representation 
matrix that would act on the 2 N dimensional space of degenerate ground states were we 
to make the braid group act by holonomy on the states, instead of by explicit monodromy 
on the Hamiltonian. The 0(2iV) spin representation is irreducible, but under restriction to 
SO(2iV) it decomposes into two irreducible components, these being the spaces of odd or 
even fermion number. For the elementary braiding operation Ti, Ivanov showed Q| that 
the relevant unitary operator U (Tj) can be taken to be 

r i = -^(l + 7 i+ i7i), (18) 

as this is unitary and obeys 

TiliT' 1 = 7i + l, 
Tili+lT' 1 = -Ji. 

12 



The operation of taking vortex i completely around vortex i + 1 and back to its starting 
point is therefore 

T~hi T i 2 = -1%, 

T hi+i T ^ 2 = -7*+i- 

For the rest of this section we will take the holonomy point of view — i.e. the Hamiltonian 
will be kept fixed as Hq = £V Eib\bi, and a braiding T will act on the state by U(T). 

We now consider the effect of various braid group generators on the occupation number 
basis states 

\n 1 ,...,n N ) = (b\) ni ...(b* N ) n »\0). (19) 

When we expand out the generator Tj in terms of the annihilation and creation operators, 
we find that 

n = ^=(1 + b i+1 bi + b\ +1 h + b i+1 b\ + b\ +1 b\), (20) 

and so 

T i \n 1 ,...,ni,n i + 1 ,...,n N ) = -j= 



{K,.. 


.,n h n i+1 , . . .,n N ) 




+ |ni, . 


. . , (rii - 1), (rii+i - 1), . 


• • , n N ) 




. . ,(rij- + 


■n N ) 


- ni, . 


. .,(rij + 1), (rii+i - 1),.. 


. . , n N ) 




..,(Ui + 1), (lli+i + 1), . 


..,n N )} 



(21) 

Here we understand that when the ±1 takes the occupation number rij out of the set {0, 1} 
the illegal state is to be replaced by zero. 

We next define to be the operation of interchanging vortex i with its sibling antivor- 
tex. (A vortex and an antivortex being distinguishable, this operation does not return the 
system to its original configuration and hence cannot be considered as a braiding operation.) 
The corresponding unitary operator is 

= ^(1 + <(&&-&,&!)). (22) 

From this, we find 

tv (0 Vi, ...,m,...,n N ) = e i7Tni/4 \n!, . . . ,n h . . . ,n N ). (23) 
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One can also consider an operation involving interchanging vortex i with the antivortex of 
its partner, i + 1, to the right. We will call this operation . The corresponding operator 
is 

rl 1] = - + i&kfc - jfc^&J + j6 t fi6 t) > (24 ) 

and this acts as 

T^\ni, . . . ,ni,n i+1 , . . . ,n N ) = — = {\m, . . . , n», . . . , n N ) 

-i\ni, ...,(rii- 1), (n i+1 - 1), ... , n N ) 
+i\n 1 , . . . , (m - 1), (n i+1 + 1), . . . n N ) 
+i\n 1 , . . . , (m + 1), (n i+1 -l),...,n N ) 
-i\nx, ...,(m + 1), {n i+ i + 1), . . . ,n N )}. (25) 

A strategy for accessing the protected information is now apparent. The wavefunction of 
the protected state |\&) belonging to the even (or odd) fermion number sector is specified 
by the 2 N ~ 1 complex numbers forming the its components in the occupation-number basis 
\ni, . . . , tin) . Provided that we can have access to multiple copies of we can repeatedly 
fuse the vortex-antivortex pairs and so estimate the probability of a particular pattern 
(ni, . . . ri/v) of relict particles. In this way we obtain the numbers |(\l/|ni, . . . ,n7v)| 2 - We lose 
all relative phase information in this process, however. This loss occurs not only because 
we are finding probabilities, but also because the different occupation-number states have 
unpredictably different tunneling energies, and so time evolution will scramble their relative 
phases during the fusion process. 

Not all is lost, however. (|21j) and ()25j) seem to suggest that we can use direct fusion 
combined with controlled braiding r< and interchange t$ to find all three of the numbers 



l^oil 2 


= |<*K. 


..,0,1,.. 


•,^7V)| 2 , 




lAol 2 


= |(*|m,. 


..,1,0,.. 


■,n N }\ 2 , 




\A 10 + A 10 \ 2 


= 


..,0,1,.. 


.,n N ) + (*|ni, . . 


.,l,0,...,n w )| 2 , 


\A 01 +iA 10 \ 2 


= |(*K,. 


..,0,1,.. 


.,n N ) + i(^\n 1 ,. 


. .,1,0, ...,n N }\ 2 



If this is possible, we would be able to obtain all relative phases of (^\rii, . . . ,njv)'s, for 
given \zi\ 2 , \z 2 \ 2 , \zi + z 2 \ 2 and \z\ + iz 2 \ 2 , it is possible to recover the complex numbers z 1 
and z 2 up to a common phase factor. 
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This is not quite the case, however. A vortex-antivortex interchange would leave one 
vortex- vortex pair and one antivortex-antivortex pair. These pairs are problematic, because 
the excited state of this pair would not lead to an unbound Bogoliubov quasiparticle exci- 
tation, but rather a bound excited state of a 'double half-quantum' vortex (or antivortex). 
This means that fusing such pairs would not give us a result that we can read out. In short, 
\z\ + 1Z2 1 2 is not available to us. 

Nevertheless, there is still much we can figure out about relative phases of 
(^\fii, . . . , 7ijv)'s. For iV complex numbers, z^ if one knows all their absolute values and 
distance between each other, \zi — Zj\, for all i,j (or equivalently \zi + Zj\), the geometric 
configuration of determined rigidly. By this we mean that we have determined z^s 

up to overall rotation around the origin and reflection with respect to the real axis (that 
is, complex conjugation). It needs to be noted that knowing \zi ± Zj\ for all i,j actually 
overdetermines z^s. So it turns out that even though direct fusion combined with controlled 
braiding only gives us a subset of {\zi ± Zj\}, it is still sufficient for us to determine 
(fy\ni, . . . , n7v)'s up to complex conjugation and overall phase. The exact procedure will be 
given in the appendix. (One needs to be note that if there is a superposition of even total 
occupation number states and odd total occupation number states, there is no way to obtain 
phase relation between them, which reflects the fact that, as Ivanov pointed out [17|, the 
superconducting Hamiltonian creates or destroys electrons only in pairs.) 

IV. CONCLUSIONS AND OPEN QUESTIONS 

Half-quantum vortices in the p x + ip y superconductors provide, at least in principle, a way 
to generate and manipulate entangled states in a topologically protected Hilbert space. We 
have shown that by physically fusing vortices with appropriate antivortices it is possible to 
reconnect the protected space to the rest of the Hilbert space and so read out the information 
encoded there. 

A number of questions remain, however: i) We have assumed that the fusion process 
is slow enough that the bound state merges adiabatically with the continuum. What 
happens if the process is too fast? Can we find a tractable model for the annihilation process 
that would allow us to determine how adiabatic it has to be in order not to lose information? 
ii) Our picture of generating and moving vortices is at present only a thought experiment. 
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The necessary half-quantum vortices have not even been detected in any real system. If 
they can be found, can we come up with some practical device (a configuration of STM 
probes, current sources, drains, gates etc) that can create, guide, and monitor the vortices? 
Since no experiment is likely to be able to detect a single quasi-particle, we would have 
arrange for a continuous operation and measure currents; iii) Perhaps the most interesting 
question is whether we take the insights developed from the relatively simple picture of 
braiding and fusion in the superconductor and extend them to the other candidate systems 
with non-abelian statistics. 
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APPENDIX A: BERRY PHASES AND BOGOLIUBOV TRANSFORMATIONS 

In this appendix we review the formal algebraic aspects of Bogoliubov transformations, 
and their implications for the computation of Berry phases for BCS states. 



1. BCS ground state 

Suppose that Hij is an N-bj-N matrix representing a one-particle Hamiltonian. When 
we include the effect of a superconducting condensate, the second-quantized Bogoliubov 
Hamiltonian becomes 



-^Bogoliubov — CLiHijdj + -Aij^Qj + -A^a^a 




A f . -H T 




(Al) 



Here a\ and a« are fermion creation and annihilation operators, the gap function Ay is a 
skew symmetric matrix, and H T denotes the transpose of the hermitian matrix H . For 
a continuum superconductor, the index "i" should be understood to incorporate both the 
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space co-ordinate x, and the spin index. A sum over i therefore implies both an integral 
over real space and a sum over spin components. 

The many-body Hamiltonian is diagonalized by means of a Bogoliubov transformation. 
To construct this transformation we begin by solving the single-particle Bogoliubov-de- 
Gennes (BdG) eigenvalue problem 

H A \ ( u n \ I u r 

= E a I . (A2) 




At -H 





E a \ ° ), (A3) 



At - 

Here u a and v a are A^-dimensional column vectors, which we take to be normalized so that 
\ua\ 2 + \v Q \ 2 = 1- If we explicitly display the column- vector index i they become matrices 
u ia and v ia . Taking the complex conjugate of (jA2|) tells us that 

H A \ ( v* 

K 

and so the BdG eigenvalues come in ± pairs. We will always take E a to be the positive 
eigenvalue. 
We now set 

a% = u ia b a + v* iQ b ] a 

a\ = v ia b a + u* a b ] a . (A4) 

The mutual orthonormality and completeness of the eigenvectors (u a ,v a ) T ensures that the 
b a , W a have the same anti-commutation relations as the a« a\. In terms of the b a b^, the 
second-quantized Hamiltonian becomes 

N ^ N . N 

-^Bogoliubov = ^2 Eab^ba — - ^ E a + - ^ . (A5) 

a=l a=l i=l 

Here the Ef ] are the eigenvalues of H. These can be of either sign. 

If all the E a are strictly positive, the new ground state is non-degenerate and is the 
unique state |0) 6 annihilated by all the b a . If we could find a unitary operator U such that 

b a = aiU* ia + a\v* a = Ua a U~ l 

b a = a\u ia + a i v ia = Ua) a U~ 1 (A6) 

then we would have |0) 6 = U\Q) a , where |0) o is the no-particle vacuum state. It is not easy 
to find a closed- form expression for U, however. An alternative strategy for obtaining |0) 6 
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begins by noting that if that the matrix U{ a is invertible then the condition 6j|0) b = is 
equivalent to 

(a, + 4 u fca(w* _1 )«)|0) 6 = 0, i = 1, . . . N. (A7) 
We therefore introduce the skew-symmetric matrix 



Sij = OX 1 (A8) 



and observe that 



exp { \a\a]Sij 1 a k exp j -^ajajSy J> = a k + a\S ik . (A9) 



2 

From this we conclude that 

|0) 6 = iVexp|iaJa]5 iJ j|0) a (A10) 

where |0) a is the original-no particle state. Equation ()A10|) explicitly displays the super- 
conducting ground state as a coherent superposition of Cooper-pair states, and allows us to 
identity with the (unnormalized) pair wavefunction. 
The normalization factor M is found from 

(S^) = det 1/2 (I + S\S 2 ), (All) 

where 

\S) = exp {-ajaj^j |0) a , (A12) 

to be 

Ar = det- 1/4 (/ + 5 t 5). (A13) 

Because the group of Bogoliubov transformations on a iy a\, i = 1, . . . , N is SO(2N), and 
because the subgroup U{N) ~ Sp(2A^, K) R SO(2A^) of transformations that mix the a's only 
with themselves (i.e. not with the a^'s), preserves the no-particle vacuum |0) a , the set of 
physically distinct ground states is parameterized by the symmetric space SO(2N)/XJ(N). 
As a check of this assertion, observe that 

dim SO(2A0 - dim U(JV) = N(2N - 1) - N 2 = N(N - 1), (A14) 

which is the number of independent real parameters in the complex skew-symmetric matrix 
S^. These for i < j serve as complex co-ordinates on all but a set of measure zero in the 
manifold of possible ground states. 
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2. Clifford algebra and Lie(S0(2iV)) 



We can make the S0(2iV) character of the Bogoliubov transformations manifest by in- 
troducing a set of 2N Dirac gamma operators. These are related to the fermion annihilation 
and creation operators by 

7i = (oi + aj) 
7 i+A r = i(aj - Oj). (A15) 

The 7j are hermitian, and obey the Clifford algebra 

H, lj} = lilj + Ijli = 2S ij- ( A16 ) 



The Hamiltonian can be rewritten in terms of the 7, as 

2N 

1 

-^Bogoliubov 



^E^. + itrF (A17) 

where 



r« = ^[7i,7i] (A18) 

are the spinor generators of the Lie algebra of SO(2AQ, and the matrix hij has entries 

/ -SfH - 3A -£if + KA \ 
/ii 7 - = . (A19) 

y + ka + ^a y 

Here = X and = Y denote the real and imaginary parts of Z = X + iY, and the 
vector space has been partitioned so that i — 1, . . . , AT is the first block and i — N + l, . . . 2N 
is the second. This rewriting reveals that a general Bogoliubov Hamiltonian is, up to an 
additive constant, an element of the Lie algebra of SO(2iV). The Bogoliubov transformation 
that diagonalizes H is therefore the operation of conjugating the Lie algebra element H into 
the Cartan sub-algebra. This we can take to be spanned by the commuting set of operators 

^ [jn+z, 7»] = 2 fata _ a * a l ) ' i = 1 ' • • • ' N - ( A2 °) 

After conjugation, 

-^Bogoliubov ~~ ¥ ^-f^Bogoliubov^ \ (A21) 



with 



u = ex p U E mv> [ ( A22 ) 



2 
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for some parameters the Hamiltonian becomes. 

N 1 1 N 

-^Bogoliubov — > ^ -Bq — [7jy +a , 7q,] + ^ ^ -E^ 
i=l i=l 

1 N 1 N 

= 2E £tt ( a « flft-M «) + 2E £,(0) 
i=i i=i 

N 1 1 

= ^^aK-^E^+oE^' ( A23 ) 



which is the same as (jA5|) . Again it is convenient to regard the energies E a as being positive, 
even though the E®, the eigenvalues of H, can have either sign. The additive constant is 
then the ground-state energy of superconducting system. 



3. Zero modes 



The gamma-operator language is particularly useful when there are zero modes. Because 
of the ±E symmetry, any zero energy eigenvectors of the BdG Hamiltonian must come in 
pairs. Suppose that (u ,v ) T becomes degenerate with its negative energy sibling (vq,Uq) t . 
Then we can write this pair of eigenvectors' contribution to the mode expansion as 



where 




(A24) 



(A25) 



and 

lo = (b + bl), lN = i (bl-b ). (A26) 

The column vectors (U , V ) T and (Un,Vn) t both have the feature that U* = V. This 
anti-linear up-down symmetry is characteristic of the localized zero modes in the vortex 
cores. 



20 



4. Berry Connection 



We need to compute the Berry connection iA = (S\d\S), where \S) = Af\S) is the 
normalized ground state. To do this we exploit the fact that the un- normalized states \S), 
being functions only of the <Sy, and not of the S*j, defines a holomorphic line-bundle over 
the Kahler manifold SO(2iV)/U(iV). We can therefore read off the 1-form connection iA 
from derivatives of the Kahler potential 



as iA = dlnAf - dlnAf 



lnjV= -~ lndet (I + SflS) (A27) 
If we express the parameters 5V,- in terms of the normalized 



Bogoliubov eigenvectors (u, v) T , we find that 

91nA/, 

rV.SV. 11 OS,;, 



_ f dlnAf dlnAT 

i<j y % 3 L J 



, N / v * \ 

= oX)^" Ua)d \ * +^{ Ar g( detM )l- (A28) 

a=l \ u a J 

The expression f|A28j) has a simple interpretation. From ()A3|1 we see that the column vec- 
tors (f*,w*) T are the negative-energy eigenstates of the one-particle Bogoliubov-de-Gennes 
Hamiltonian 

U -%)■ <a29) 

If this were a Dirac fermion problem, we would fill the Dirac sea consisting of these negative 
energy states. The Berry phase of the vacuum would then be the sum of the Berry phases 
of the occupied states. The first term in (|A28|) is precisely one-half of this sum. The factor 
of one-half compensates for the artificial doubling of the degrees of freedom in passing from 
H to i^BdG- The second term in (jA28j) is a total derivative, and reflects a choice of gauge. 
We next compute the Berry connection for an excited state 

\a 1 ,...,a n ) = bi i ---bijS), (A30) 

by using 

dbl = (ui/sbfs + v*pVp) dv ia + {v ip b p + u* if3 b^) du ia 

= {v*pdv ia + u* if3 du ia ) tfp + (Uipdv ia + Vipdu ia ) bp. (A31) 



21 



When the state |«i, . . . , ot n ) is non-degenerate, we are interested only in the diagonal (3 = a 
term 

c^tjdiag = (v* a dv ia + u* a du ia ) bl, (no sum on a), (A32) 

and we find 

iA = (ati, . . . ,a n \d\cti, ...,a n ) 
= (S\b an ---b ai d{bl---bl \S)) 
= (S\d\S) + (S\b an ---b ai d(bl---bl n )\S) 

= (S\d\S) + J2« m <Jd( Uam ), (A33) 

m=l \ V am J 

which is the sum of the many-body ground-state Berry connection and the Berry connections 
of the individual one-particle excited states. Observe that there is no factor of "1/2" in the 
contribution of these occupied excited states. 

The non-abelian Berry connection of a set of degenerate many -body states is computed 
in the same manner. It will include a diagonal term from the reference state \S) and a sum 
of non-diagonal terms terms of the form 

f u a \ 

iA Pa = {u* p ]■ ( A34 ) 

In the case of the vortices, the states of interest are the exponentially localized core states. 
Because of this localization only the overlap of each core state with itself has any chance of 
providing a non-zero term in the connection, but it is readily verified that with the phase 
choices made in the text, all these contributions are zero. The only "Berry phase" produced 
ry the vortex braiding is the overall diagonal Berry phase associated with the Magnus force 

APPENDIX B: DETERMINING PHASE RELATION AMONG COEFFICIENTS 
OF OCCUPATION NUMBER BASIS 

Let us consider how many vortex interchange steps would be needed in order to rigidly 
- this rigidity being defined in the last page of the section III - determine geometric con- 
figuration on the complex plane of all {^\ni, . . . , n^'s, with having a definite parity in 
total occupation number. 
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In the case N = 2 one can easily see from 

r 1 (A oo |00)+A 11 |ll)) = -L{(A 00 + A 11 )\00)-(A 00 -A 11 )\11)} : 

n(A 01 |01) + Aio|10» = ^={(An + Ai )|01> - {A 01 - A 10 )|10>} (Bl) 

that the fusion following the vortex interchange process gives us, up to the sign, phase differ- 
ence between A 00 and A u (and likewise between A Q1 and A w ). Given that we already know 
\A Q0 \ and \ Au\, this is sufficient to rigidly determine on the complex plane the configuration 
of Aqo and An. Same can be said for A$i and Aiq. 

It is instructive to work out the next simplest case N = 3, where |\&) = A oo|000) + 
Ano 1 110) + Am|011) + Aioi|101). From N = 2 case, one can see that T\ rigidly determines 
the configuration of A oo and Au on the complex plane. It also rigidly determines the 
configuration of A ou and A 101 on the complex plane. However any phase relation between 
the former and latter remains completely unknown at this point. But then from 

r 2 |tt) = -4{(^ooo+Aoii)|000)-(A 1 oi-A 110 )|110)-(Aooo-Aoii)|011) + (Aioi+A 1 io)|101)}, 
V2 

(B2) 

one can see that after fusion following the implementation of r 2 , the configuration of both 
the A 00Q , A u pair and the A W i, A\ W pair would be determined rigidly with respect to the 
origin. We now have a rigid configuration of all four coefficients - Aooo? Aon, Aloij and Auq 
- on the complex plane; they have been determined up to complex conjugation and overall 
phase. 

The N = 3 case gives us ideas about how to make use of a recursion argument for the 
general case in figuring out the phase relation. One can assume that for iV = m there is 
some process consisting of fusion of vortex- antivortex pairs and applications of r^'s (where 
k < m — 1) which give us the rigid configuration on the complex plane of the coefficient of 
states with even total occupation number. The same process would also give us the rigid 
configuration on the complex plane of the coefficient of states with odd total occupation 
number as well. (We have seen that this holds true for m = 2.) Now note that the even (or 
odd) total occupation number sector of the Hilbert space in N = m + 1 case can be divided 
into the following two classes: 

|ni, . . . ,n m ,0), (B3) 
K,...,rCl>. (B4) 
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The operations that were used in the procedure we have applied for obtaining all phase 
relation in the N = m case does not affect n m+ i, and so by applying these operations we 
would obtain the rigid configuration for the coefficients of states belonging to (jB3|) . Same 
can be said for the coefficients of the states of belonging to (jB4J) (though |ni, . . . ,n m ) and 
\n[, . . . , n' m ) have opposite parity in total occupation number). Now all we need to do is to 
figure out the rigid configuration of two pairs of coefficients, each of which consists of one 
coefficient for one of the states belonging to (jB3|) and one coefficient for one of the states 
belonging to (jB4j). So for |*) = J] An,...,n m ,n m+ i • • ■ ,n m ,n m+1 ), (note that summation 
is restricted to even, or odd, total occupation number) 

= ~/| { ( Ati , „ .,n m _ i ,0,0 + ^m,...,n m _i,l,l)Kl, • • • ,™ m -l,0,0) 

— (Aii,...,n m _i,0,0 — Aii,...,n m _i,l,l)|^l, • • • , 1; 1)} 
+ ^(^^Mr-.n'm-lfiA + ^n' 1 ,..,n^_ 1 ,l,o)|^'l, • • • , 0, 1) 

- (^,...,^,0,1 - A n ' 1 ,..,n^_ x ,i,o)|ni, • • • , 1, 0)}. (B5) 

One can easily see that all pairs belonging to (An,...,r»m-i,o,o» An,...,n w -i,i,i) or (A»i,...,n^_ 1 ,o,i, 
^4ni,...,n' _ 1 ,i,o) now have rigid configuration. It is clear that we now have rigid configuration 
for {A 

m,...,n m _i,n m ,n m+ i}- Also we can see that in this scheme of figuring out the rigid 
configuration of 2 N ~ 1 coefficients of even (or odd) fermion sector in the occupation-number 
basis, N — 1 vortex interchange steps are needed. 
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